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Abstract: The influence of impurities, embedded into the isolating spacer (I) 
between two ferromagnetic electrodes (F), on the I-V curve and tunnel magnetore- 
sistance (TMR), is theoretically investigated. It is shown, that the current and TMR 
are strongly enhanced in the vicinity of the impurity under the condition that the 
energy of the electron's bound state on the impurity is close to the Fermi energy. If 
the position of the impurity inside the barrier is asymmetric, e.g. closer to the one 
of the interfaces F/I the I-V curve exhibits quasidiode behavior. 

The magnetic tunnel junction (TMJ), consisting of two metallic ferromagnetic electrodes 
separated by insulating barrier and exhibiting tunnel magnetoresistance (TMR) of order 50% 
attracts a lot of attention □ □ □ especially due to their_ possible application in MRAM 



he theory of TMR for the ideal 
that in presence of different 



(Magnetic Random Access Memory). In the pioneer paper 
(without defect) TMJ was developed. Later it was shown [5, 
types of defects within the barrier the I-V current and TMR change dramatically. In these 
papers the averaged over cross section of the system current was calculated. However it is 
interesting to investigate the local (in vicinity of the impurity) current and TMR, especially 
taking into account that using the STM technique it is possible and was realized [tJ to 
span tunnelling current over the cross section of TMJ. Recently the theory of local impurity 
assisted tunnelling in TMJ was developed [8]. As a model of TMJ was taken tight binding 
model and Kubo formalism was used for calculation of spin-dependent tunnel current. In 
this paper the I-V curve was not investigated in detail and besides that the dependence 
of spin-dependent current on the position of cross section plane relative to the position of 
impurity was not investigated. 

In the presented paper we investigate the local distribution of spin-dependent current 
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for different positions of the cross section plane and local I-V curves of TMJ with single 
impurity and for random distribution of impurities inside the barrier. We adopted the free 
electron model with exchange splitting for ferromagnetic electrodes and used nonequilibrium 
Keldysh technique ji| for calculating of nonlinear on applied voltage transport properties. 

We considered the model of TMJ as a three layers system, consisted from two thick 
ferromagnetic electrodes F separated by the insulating layer, I . Inside the barrier the single 
nonmagnetic impurity with attracting potential was situated at some distance from F/I 
interface. The two cases were investigated : parallel and antiparallel orientations of F-layers 
magnetization. 

The F-electrodes are connected to the reservoirs with chemical potentials yUi and \i 2 so 
that fj,2 — A*i = eV, where V is the applied voltage. 

To calculate the current through the system we have to found Keldysh Green function 
G h and advanced and retarded Green functions G A and G R . Solving the Dyson equation 
we found that 



where Gq + (t, r'), G A (r, r') and G R (r,r') are the Green's functions for the system in the 
absence of the impurity and the potential of the impurity V was represented as ^-function: 
V(r) = Wal5(z — z )5(p-p ), r = (p , z ) is the position of the impurity, a is it's effective 
radius, W is it's intensity. The explicit expressions for G A , G R , G~ + have the following form: 



G- + (r,r') = G + (r,r / ) + 



Gff(r,r )iyG - + (r ,r') G„+(r, r )WG$(r , rQ 
l-^(r 0l r ) l-WG A (v ,r ) 

G R (r, r )^C - + (ro, r )^(r , r') 
+ (1 - WG R (r , t )) (1 - WG A (r , r )) 



(1) 




{z'iZx) [q(zi) - ih)} , 



{E(z 2l z) [q(z 2 )+ik 2 } + E- 



(z 2 ,z) [q{z 2 ) - ik 2 ]} 



(2) 




{E{z 2 ,z) [q(z 2 ) - ik 2 ] + E- 1 



(z 2 ,z) [q(z 2 ) +ik 2 ]} 



(3) 



x {E{z',zx) [q{z\) — iki] + E~ 



{z',zi) [q(zi) +ih]} 
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G +(r,r') = j^ R 1^f^B±l {E(z',z 2 ) [q(z 2 ) + ik 2 ] + E-\z> , z 2 ) \q{z 2 ) - ik 2 ]} 

W q(z)q(z')\den\^ 

x {E(z, z 2 ) [q{z 2 ) - ik 2 ] + E~\z, z 2 ) [q{z 2 ) + ik 2 }} 

(4) 

x {E(z u z) [q(z x ) - ih] + E~ x {z x ,z) \q{z x ) + ih]} , 
where 



q(z) = x /qi + ^-^{^- ) eV, 



k^J^ie-A,)-^ 



den = {E(z 1 . i z 2 ) [q(z 2 ) - ik 2 ] [q{z x ) - ik x ] - E~ 1 {z l , z 2 ) \q{z 2 ) + ik 2 ] [g(zi) + ih]} , 
E( Zl ,z 2 )=e^ q{T)dT , 

n is the electron momentum perpendicular to the plane of structure, e is the energy, z\ and 
z 2 are the positions of F/I interfaces, Ai and Ai denote the positions of energy band bottom 
for spin up and down subbands. 

til = and tir = f°(e + eV) are Fermi distribution functions in the left and right 

reservoirs and -g- 2 - height of potential barrier above Fermi's level. 



In (01,(121),© and ® p and z are in the plane and perpendicular to the plane coordinates, 
and we consider that z and z are situated within the barrier. We have to take into account 

t I t I 

that all Green functions are matrixes in spin space. We have consider k[ F , k\ F , k 2F) k 2F are 
Fermi wave vectors of electron with spinf (j) in the left and right F-electrodes. The current 
was calculated, using the following expression: 

eft f ( dG~+(z, p; z',p) dG- + (z,p; 
3 ^ Z) -^J d£ { dP dz ) z=z , (5) 

On the FigQ and Fig|21 the dependencies of the currents in different channels ( up and 
down spin) on coordinate p — po at one interface I/F (z 2 = 15 A) (another interface is at 
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FIG. 1: Dependence of the current for different spin channels and P and AP configuration on the 
distance from the impurity in the plane of the structure at z=15 A. k F = 1.1 A , kp = 0.6 A , 
go = 1.0 A" 1 



z\ = 0) and inside the barrier at z = 10 are shown. Position of the impurity is at p — 
and z = 5 A. 

So in the vicinity of the impurity the hot spot of the radius approximately equal 6A may 
be observed and the value of the current density in the center of the hot spot exceeds the 
value of the background current on several orders of magnitude. On the FigEl the TMR 
dependence on the distance from the impurity at different z is shown. It is interesting that 
the value of TMR in vicinity of the impurity exceeds it's background value (TMR for the 
ideal structure is equal 0.013) more over then order of magnitude and for some cases it exists 
the region of p — po where TMR becomes negative. 

Now on FigEJthe I-V current for positive and negative applied voltage is shown. This 
curves are quite asymmetric on the sign of the voltage. It is connected with asymmetry of 
position of the impurity inside of the barrier. We choose the potential of impurity so that 
bound (resonance) state of electrons with spin up located near Fermi level for the positive 
applied voltage = 1.2 V, and for negative voltage the position of this bound state lies below 
Fermi level. This diode behavior is demonstrated for the single impurity and to have the 
possibility to use this diode effect in practice we have to investigate the case of the finite 
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FIG. 2: The same dependence at z = 10 A. 
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FIG. 3: Dependence of TMR on the distance from the impurity in the plane of the structure at 
different z. For parameters see FigQ 



concentration of impurities. 

In this case we consider the same magnetic tunnel barrier structure with monolayer of 
impurities of finite atomic concentration x, situated closer to the one of F/I interface. To 
solve the problem as a first step we have to find coherent potential and effective Keldysh 
Green function G~g . Solving the Dyson equation in the Keldysh space we got the following 



6 





50 -i 




45 - 




40- 




■ 

35 - 




30- 






cm 


25- 


$ 






20- 




15- 




10- 




5- 




0- 



Left to Right 
Right to Left 



o.o 




V, Volt 



FIG. 4: Local I-V curve at p = po and z = 15 A for the case of single impurity. 
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FIG. 5: I-V curve in the case of the layer of impurities at zq = 3A and x = 0.5. 



expression for G+ + ' ' ' 



ft 



G~ + (z z) = G + (z z) + ^° + [h ^ AG o (fpj z ') _|_ ( z i z o) ^ Rg o + ( z o, z') 



l-G$ (zq, zq) T, a 1-G*(zo,zo)E r 

Go(z,Zo)T, + Gq(zq,z') Gq (z, Zq) S^Gg + (Zq, Zq) S A Gq (Zq, z') 



(6) 



(1 - G$ (z , z ) E A ) (1 - Gfi (z , z ) £«) (1 - Gfi (zo, z ) Z A ) (1 - Gq (z , z ) £*) 
where are the coherent potential (CP.) for the retarded and advanced Green functions, 
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which has to be found from the C.P.A equation: 

t=(l-x) ^ - ^ + (x) ^ - ^ = (7) 

1 - {e A - X)G eff (z ,Po; z ,p ) 1 - (e B - E)G eS (z , p ; z Q , p ) 

where e A and e B are the onsite energies of the host (AI2O3) and the impurity (Al) and 
Z-+ = i(n R + n L )(£ R -X A ). 

Now to calculate I-V curve we may use the found G aa , substituting it into the 
expression flJJJ). 

On the FigJ5] the I-V curve for AP configuration is shown, and the asymmetry of the 
curve on the sign of applied voltage is clearly seen. 

Such a structure may be prepared if to sputter thin layer of Al on the F-electrode, then 
oxidise it, after sputter thicker layer of Al and oxidise it from the top side not completely. 
So some thin layer of the random alloy Al x Al20$(\- x ) is situated inside the more or less ideal 
insulator Al 2 0^ at the distance close to the first F/I interface. 

The work was partly supported by the Russian fund of fundamental research (grant N 
04-02-16688a). 
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